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Abstract. We review the properties of Galactic binaries containing two compact objects, as derived by means 
of population synthesis. Using this information we calculate the gravitational wave signal of these binaries. At 
frequencies below f ^ 2 mHz the double white dwarf population forms an unresolved background for the low- 
frequency gravitational wave detector LISA. Above this limit some few thousand double white dwarfs and few 
tens of binaries containing neutron stars will be resolved. Of the resolved double white dwarfs ~ 500 have a 
total mass above the Chandrasekhar limit. About ~ 95 of these have a measurable frequency change allowing a 
determination of their chirp mass. We discuss the properties of the resolved systems. 
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1. Introduction 

The interest in gravitational waves, predicted by 
Einstein's theory of general relativity, was greatly en- 
hanced by the signals supposedly detected by resonant 
gravitational wave (GW) antennas (Weber 196£) and the 
discovery of the pulsar B1916+13 in a relativistic bi- 
nary ( |Hulse fc Taylor 1975| ; [Taylor fc Weisberg 198^ ). 
Currently, about ten projects for ground and space-based 
gravitational wave detectors are already operating or un- 



good (see Kalogera et al. 200(]| ). They could be better for 
black hole - black hole or for black hole - neutron star 



mergers (Tutukov & Yungclson 1993b; Lipunov et al. 1997 



Portegies Zwart fc McMillan 200C) 



der development (see Flanagan 199S). They will open the 
windows in the frequency bands 10 to 10^ Hz from the 
ground and 10^*^ to 1 Hz from space , Recently the first 



At low frequencies, it was first expected that contact 
W UMa binaries will dominate the gravitational wave 
spectrum (Mironovskii 1965). However, it was shown that 
the gravitational wave background formed by Galactic 
disk systems is probably totally dominated by detached 
double white dwarfs and that their number is so large that 
they will form a confusion limited background for the cur- 
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upper Limits on detections from the Japanese TAMA300 



rently planned detectors (Evans et al. 1987; Lipunov et al 
1987| ; |Hils et al. 1990| ; [Nelcmans et al. 2000aD . Only sources 



detector were reported (Tagoshi et al. 2001) 
At i high froquoncic i 



the merging ovonto of e3ctragalo 



with a frequency above a certain limiting frequency (some- 
where between ^ 1 — lOmHz) can be resolved (Evans et al 



tic bin hrico containing neutron otaro and/or black holoo 
are among the most promising sources of GW radiation. 
The estimates of the merger rates of these systems are 



1987) 



highly uncertain 


(c-g. 


Phinney 1991 




Portegies Zwart & 


Yungelson 1998 




Kalogera & Lorimer 2000|). An upper 



limit for the rate of neutron star 
in our Galaxy of 



neutron star mergers 



10 ^yr ^ is found both from obser- 



vations ( Arzoumanian et al. 1999) and theory (Tutukov 
fe Yimgelson 1993b). Extrapolated to cosmic scales, these 
estimates show that the perspectives for detection of such 
events by the first generation GW detectors are not very 
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The aim of the present paper is an accurate evaluation 
of the confusion limit, based on population synthesis mod- 
els for compact stars in the Galactic disk and a discussion 
of the properties of the sample of potentially resolved bi- 
naries containing two compact objects: white dwarfs, neu- 
tron stars or black holes. We first discuss the gravitational 
wave signal from (eccentric) binaries (Sect. |). Next, we 
summarise the properties of the Galactic disk populations 
of compact binaries which are relevant to the emission of 
gravitational waves (Sect. I). We do not consider globular 
cluster binaries. In Sect. ^ we present a model for the back- 
ground formed by the Galactic disk double white dwarfs, 
discuss the confusion limit and the properties of the in- 
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dividually resolved binaries. A discussion of the possible 
contribution of the halo and extragalactic sources and a 
comparison with previous work follows in Sect. || Our 
conclusions are summarised in Sect, pi 



2. Gravitational waves from binaries 

The gravitational wave luminosity of a binary in the nth 
harmonic is given by (Peters & Matthews 1963) 
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L{n, e) 



32 hp {M + m) 



g{n,e) 



(1) 




Here M and m are the masses of the components, a is 
their orbital separation and e is the eccentricity of the 
orbit. The function 17(71, e) is the Fourier decomposition of 
the GW signal. 

The measurable signal for gravitational wave detectors 
is the amplitude of the wave - /i+ and hx for the two 
polarisations. These can be comput ed from the GW flux 
at the |bartii (Vrcss & morne l\)('z) 



1 2 3 4 5 6 7 



9 10 11 12 13 14 15 16 17 



Fig. 1. Scale factor of the GW strain amplitude 
\/ g{n, e)/n for the different harmonics [Eq. (^)] for e = 0, 
0.2, 0.5 and 0.7. 



1996; Portegies Zwart fc Yungelson 1998; Nelemans et al 



2001b). The basic assumptions used in this paper can be 
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Assuming the waves to be sinusoidal and defining the so 



called strain amplitude as h 
obtains 
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summarised as follows. The initial primary masses are dis- 

(2) tributed according to a power law IMF with index —2.5, 
the initial mass ratio distribution is taken flat, the initial 
semi major axis distribution flat in log a up to a = IO^Rq, 
and the eccentricities follow Pie) oc 2e. The fraction of bi- 
naries in the initial population of main-sequence stars is 
50% (2/3 of all stars are in binaries). A difference with 

(3) other studies of the populations of close binaries is that 
the mass transfer from a giant to a main sequence star of 

-1 comparable mass is calculated using an angular momen- 
7 tum balance formalism, as described in Nelemans ct al 



(2000t). For the star formation rate of the Galactic disc 



where M. = (Mrtif^^ l{M -I- m)^/^ is the so called chirp 
mass and cjg ~ 'nnjPoY\i is the angular frequency of the 
emitted wave^. In Fig. || we plot the values of ^ g{n, e)/n 
for different eccentricities. High eccentricity binaries emit 
most of their energy at higher frequencies than their or- 
bital frequency, reflecting the fact that the radiation is 
more effective near periastron of the orbit. Thus, eccen- 
tric compact binaries may be detectable sources of GW 
signals at frequencies higher than their orbital frequency 
(cf. [Barone et al. 1988|; |Hils 199lD. 



we use an exponential function: 
SFR(t) = 15 exp(-VT) Mq yr" 



(4) 



where r = 7 Gyr. With an age of the Galactic disk 
of 10 Gyr it gives a current star formation rate of 
3.6 M(7) yr~"^ compatible with observationa l estimates 
(Rana 1991; van den Hoek fc de Jong 1997). It gives a 



3. Tlie Galactic disk population of binaries 
containing two compact objects 

We calculated the Galactic disk population of bina- 
ries containing two compact objects using the popula- 
tion synthesis code SeBa (Portegies Zwart fc Verbunt 



For circular orb its this equation is identical to Eq. (5) of 



Evans et al. (1987 ). It is different by a factor of -\/8 from 



Eq. (20) of press fc Thorne (1972| ), who use a factor \/2 
larger definition of h and possibly confuse ujg in Eq. (kI) with 
the orbital a ngular frequency. It differs by a factor 2^" from 



Galactic supernova Il/Ib rate of 0.02 yr^^ and if super- 
novae la are produced by merging double carbon-oxygen 
(CO) white dwarfs it gives a Galactic rate of 0.002 yr~^. 
Both are in agreement with observational estimates by 
[Cappcllaro et al. (1999| ). 

The current birth- and merger rates and total num- 
ber of systems in the Galactic disk with these assump- 
tio ns are given in Table p]. We use a notation introduced 
by Portegies Zwart fc Verbunt (1996| ): wd, ns and bh for 
white dwarf, neutron star and black hole respectively; ( ) 
and [ ) for detached and semi-detached binaries. The 
fact that the numbers here are different from the num- 
bers given in Portegies Zwart fc Yungelson (1998| , |1999| )p| 



Eq. (3.14) of Douglas fc Braginsky (1979 ) because they confuse 
the orbital frequency in their Eq. (3.13), with the frequency of 
the wave (twice the orbital frequency) in their Eq. (3.4). 



Note that Portegies Zwart fc Yungelson (199S| ) consider 
only a subset of the (ns, wd) population, namely the systems 
that are eccentric and contain a white dwarf more massive 
than 1.1 Mq. For this subset they find a birth rate that is 
comparable to the birth rate of double neutron stars. 
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Table 1. Current birth rates (t^) and merger rates (j^mcrg) 
per year for Galactic disk binaries containing two compact 
objects and their total number {^) in the Galactic disk, 
as calculated with the SeBa population synthesis code (see 
text). 
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2.8 X 10*^ 



Nelemans et al. (200 IbD and Brown et al. (2001) is caused 
by the differences in the assumed IMF, initial binary frac- 
tion and star formation history. 

In Fig. H we show the period distributions of the bi- 
naries of different types in the range of interest for space 
based gravitational wave detectors. The properties of these 
populations can be summarised as follows: 
Detached double white dwarf binaries: (wd, wd). 
Our model for the Galactic disk population of double 
white dwarfs is described in detail in Nelemans et ; 
(2GGlb). Most double white dwarfs have a mass ratio 
around unity and low-mass (M < 0.45Mq) components. 
From Table ^ and Fig. ^ it is clear that they vastly out- 
number all other binaries with compact objects in the 
Galactic disk. 

Semi-detached double white dwarfs (AM CVn 
stars): [wd, wd). We include in our calculation both AM 
CVn stars descending from detached close double white 
dwarfs and from low-mass helium stars with white dwarf 



companions ( Nelemans et al. 2001a ). We use Model II of 
Nelemans et al., which is most favourable for the forma- 
tion of AM CVn's. 

Neutron star white dwarf binaries: (ns, wd). 

Neutron star - white dwarf binaries fall into two fam- 



ilies ( Tutukov fc Yungelson 1993a Portegies Zwart fc 



Yungel son 1999| ; |Tauris fc Sennels 2000 ) . In one family the 
neutron star is formed first. Later the secondary forms a 
white dwarf and in the mass transfer event the orbit cir- 



cularizes (e.g. van den Heuvel & Taam 1984). If both com- 
ponents of the initial binary are of comparable mass it can 
happen that the primary becomes a white dwarf, while the 
secondary accretes so much mass that it becomes a neu- 
tron star (e.g. Tutukov fc Yungelson 1993a| ). In this case 
the orbits are eccentric. The masses of the white dwarfs 
are typically low in the first family and high in the second 
(see Fig. ^ below). 

Double neutron star binaries: (ns, ns). The for- 
mation and characteristics of the current population of 
double neutron stars is extensive ly studied by us in 
Portegies Zwart fc Yungelson (199*^ ) . Maybe the most im- 
portant assumption, which influences the birth rate, or- 
bital periods and eccentricities of neutron star - neutron 
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Fig. 2. Period distribution of the binaries of different 
types in the period range of interest to the space-based 
gravitational wave detectors like LISA. The vertical dot- 
ted lines give the periods at which the frequency of the 
fundamental (n = 2) harmonic of the gravitational wave 
is 1 and 0.1 mHz respectively. 

star binaries, is the kick velocity distribution. We use the 



one proposed by Hartman (1997). 

Black hole binaries: (bh, wd), (bh, ns) and 
(bh, bh). The knowledge of the way in which black holes 
are formed and the range of masses of their progenitors 
are still highly uncertain (see e.g. Woosley & Weaver 1995; 



Portegies Zwart et al. 1997; Ergma fc van den Heuvel 1998 



Wellstein fc Langer 1999; Fryer 1999). The treatment of 



the formation of black holes implemented in the present 
study is described in some detail in the Appendix. Typical 
black holes in our model have masses between 5 and 7 Mq. 
In the short orbital period range (Fig. ||) they are rare and 
their merger rates are at least an order of magnitude lower 
than those of the neutron star binaries (Table |]). Double 
black hole binaries are absent in this period range and do 
not merge at all in our model. 

4. The gravitational wave signal from compact 
binaries in the Galactic disk 

Merging of binaries containing neutron stars and black 
holes in distant galaxies could give measurable signals in 
the high frequency detectors. We do not extrapolate our 
merger rates to extragalactic scales, but our inferred rates 
(Table |l|) are consistent with the (very uncertain) rates 
derived elsewhere for the G alaxy (see for a detailed dis- 
cussion [Kalogera et al. 2000 ) . 

The Galactic binaries with periods less than 10 hr are 
interesting for the low-frequency GW detectors. We cal- 
culate the expected signal for LISA, the joint ESA, NASA 
detector that is expected to be launched around 2010. It 
will consist of three satellites, 5 million kilometres apart, 
between which laser beams will be exchanged, measuring 
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Fig. 3. Left: GWR background produced by double white dwarfs (both detached and semi-detached). The assumed 
integration time is 1 yr. The 'noisy' black line gives the total power spectrum, the white line the average. The dashed 
lines show the expected LISA sensitivity for a S/N of 1 and 5. Right: The number of systems per bin on a logarithmic 
scale. The contribution of the semi-detached double white dwarfs between log / ~ —3.4 and —3.0 is clearly visible. 



the distance changes (McNamara et al. 200C). It will give 
the GW amplitude as a function of frequency with fixed 
frequency resolution. A limited angular resolution will be 
achieved, allowing e.g. identification of sources in globular 



clusters (see Benacquista et al. 2001). In our calculations 
we restrict ourselves to the sensitivity in frequency and do 
not consider the angular resolution. 

Because the number of Galactic binaries drops strongly 
towards shorter periods (Fig. ||) the number of sources per 
frequency bin for detectors with a fixed frequency resolu- 
tion will also decrease: at low frequencies the signals in 
particular frequency bins will overlap, forming a so called 
"confusion limited noise" . Above a certain limiting fre- 
quency, called the "confusion limit" , there is not more 
than one system per frequency bin, so the systems can be 
resolved individually. We discuss these regimes separately. 



the resulting confusion limited background signal and the 
number of systems per bin. The contribution of the semi- 
detached double white dwarfs, which are less numerous 
than the detached double white dwarfs and have lower 
strain amplitude is concentrated in a relatively small fre- 
quency interval between log / ~ —3.4 and —3.0 where they 
dominate the number of systems per bin. 

Most previous studies used only the average properties 
of the double white dwarf population to calculate the av- 
erage background signal. In Fig. ^ we plotted the average 
of our model power spectrum as the white line. Note that 
in many bins the actual power is much larger than the 
average. These bins contain one system that has a much 
stronger signal than the rest, for example because it is 
close to the Earth, and it may be detectable above the 



noise (see Sect. 4.3) 



4.1. The confusion limited background due to double 
white dwarfs 



Evans et al. (1987) have shown that for space-born de- 
tectors the confusion limit is determined by the Galactic 
close binary white dwarfs. In our model the total number 
of detached and semi-detached double white dwarfs in the 
Galactic disk is 1.5 x 10^ (see Table |^). We distribute these 
systems randomly in the Galactic disk according to 



p{R, z) = po er^l^ sech(z//3)2 



pc 



(5) 



Here H ^ 2.5 kpc ( ^ackett 19971) and (3 = 200 pc, ne- 
glecting the age and mass dependence of (5. All systems are 
circular and for each system we calculate the strain ampli- 
tude from Eq. (||) taking Rq = 8.5 kpc and zq = —30 pc. 

To simulate the power spectrum for this population of 
binaries as would be detected by a gravitational wave de- 
tector in space we determine the distribution of the sys- 
tems over A/ = l/T wide bins, with T the total inte- 
gration time (for which we use 1 yr). In Fig. |^ we plot 



4.2. The population of resolved binaries 

Given the fact that the double white dwarf background 
buries all underlying signals at frequencies below log / ~ 
—2.8, we did not consider the neutron star and black hole 
binaries below 1 mHz. To find the binaries that will be 
resolved by LISA, we calculated the Galactic disk pop- 
ulation of all binaries containing compact objects which 
contribute to the GW signal at frequencies above 1 mHz. 
Because we now also consider eccentric binaries, which 
emit at frequencies higher than twice the orbital frequency 
[Eq. (H)], there are contributions from binaries with orbital 
periods up to '^10 hr. 

If one considers an average background as done i.e. by 



Evans et al. (1987), the average number of systems per bin 



at a certain moment drops below one system. However, we 
model all individual systems in the Galaxy and determine, 
for an integration time T of 1 yr, for each frequency bin 
(A/ = l/T) how many systems it contains. In Fig. ^ we 
plot the fraction of bins that contain exactly one, none and 
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Fig. 5. Distribution of the resolved binaries over frequency and chirp mass Ai. The grey shade gives the number of 
systems in each bin relative to the maximum in each plot which is 1735 for (wd, wd), 5 for (ns, wd) and 3 for (ns, 
ns). The type of binary is denoted at the top. The asterisk indicates the (bh, wd) system. The dots in the (wd, wd) 
panel indicate the systems with a total mass above the Chandrasekhar limit. The solid line shows the chirp line for 
T = 1 yr (Sect. |4.2| ). The curved dotted line shows the position at which the systems merge. Double neutron stars 
merge at high (/ Ri kHz) frequencies so their merger line falls off this plot. 
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Fig. 4. Fraction of bins that contain exactly one system 
(solid line), that are empty (dashed line) and contain more 
than one system (dotted line) as function of the frequency 
of the signals. For all frequency intervals the result is nor- 
malised, because the total number of bins in a logarithmic 
interval changes strongly. 



more than one system as function of frequency. The figure 
shows that the notion of a "confusion limit" as a unique 
value is too simple. At log/ = —2.84 the first resolved 
bins (i.e. containing exactly one system) are found, while 
up to log / ~ —2.3 bins containing more than one system 
are still present. 



Table 2. The number of resolved systems of different 
types (see Sect. 4.2) and the number of strong signal sys- 
tems potentially detectable above the double white dwarf 
background (see Sect. 4.3). 



Type 


resolved systems 


detectable above noise 


(wd, wd) 


12124 


5943 


(ns, wd) 


38 


124 


(ns, ns) 


8 


31 


(bh, wd) 


1 


3 


(bh, ns) 





3 


total 


12171 


6104 



The total number of resolved systems with a signal 
above the sensitivity limit (S/N = 1) of LISA for an as- 
sumed integration time of 1 yr is 12171. In Table || we 
give the numbers of binaries of different types that are re- 
solved. The eccentric binaries can contribute to more than 
one frequency [Eq. (^] , but the amplitude of the signal in 
the high harmonics in general is rather low. Just below 
the LISA sensitivity limit for T 1 yr there are indeed 3 
high harmonics of (ns, wd) binaries and one of a double 
neutron star in our model. 

The change of the frequency of a binary evolving under 
the influence of GWR is given by (e.g. Schutz 1996) 



/ = 5.8 X 1O"'^(X/M0)5/V"/^ Hzs" 



(6) 



For high-frequency systems this means that during a suf- 
ficiently long integration time the change of the frequency 
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Fig. 6. Distribution of resolved systems over the frequency and strain amplitude for the different types of binaries 
(indicated in the top right corner of each panel). The grey shade gives the number of systems relative to the maximum 
in each plot, which is 1314 for (wd, wd), 3 for (ns, wd) and 2 for (ns, ns). The asterisk indicates the (bh, wd) system. 
The dashed lines give the LISA sensitivity for an integration time of 1 yr and a signal to noise ratio of 5 (top line) 
and 1 (bottom line). The averaged double white dwarfs background is plotted as the solid line. 



(the "chirp") can be detected and hence, A4 be deter- 
mined. Then from Eq. (^ the distance to these sources 
can be found. In Fig. ^ we plot the frequency versus chirp 
mass distributions of the resolved systems for the binaries 
of different types to show their properties separately. In 
the panel with the double white dwarf systems we plot the 
so called chirp line, for T = 1 yr, as the solid line. Systems 
to the right of this line change their frequency during the 
integration by one or more bins {f T > A/ = 1/r). The 
position at which the systems merge is plotted as the dot- 
ted line in the left two panels. 

The horizontal straight dotted line in Fig. || (left panel) 
marks the lower limit of M for systems with a total mass 
larger than the Chandrasekhar mass which may be type 
la Supernovae (SNe la) precursors. The resolved systems 
above this line are plotted as the dots, since because 
of their relatively small number (501 systems) the grey 
shades above this line are practically invisible in the plot. 

In the (ns, wd) panel we also plot the chirp and merger 
line, assuming for the latter that the systems emit at the 
fundamental (n = 2) frequency. In our model the mass of 
a neutron star is between 1.25 and 1.55Mq, depending on 
the initial mass of the progenitor. This results in a very 
narrow range in chirp masses for the double neutron star 
systems. 

The quantities measured by detectors like LISA are 
the frequency and the strain amplitude [Eq. In Fig. || 
we plot the distributions of the expected resolved systems 
over log / and log h. In Fig. || we also show the sensitiv- 
ity limits of LISA for monochromatic sources, for signal 



to noise ratios of 5 and 1 and an integration time of 1 yr 
(adapted from Fig. 5 and Eq. (53) of p^jarson et al. 2000 ). 
The solid line gives the averaged noise background as pro- 
duced by double white dwarfs (see Fig. ||). 

4.3. Other detectable systems? 

It may well be that the systems which produce strong sig- 
nals because of their proximity to the Sun or a large chirp 
mass can be detected individually above the noise back- 
ground (see Fig. ||). To investigate this possibility, we com- 
puted the number of systems of the different types that 
are not resolved, but have a strain amplitude well above 
the noise background, by selecting all systems that have 
a signal above the S/N = 5 sensitivity limit of LISA (see 
Fig. H). This adds a considerable number of potentially 
detectable systems to the resolved binaries (Table §). It 
brings the total number of potentially detectable binaries 
containing neutron stars to almost 200 for an integration 
time of 1 yr. 

5. Discussion 

5.1. Halo and extra-galactic sources 

Above, the confusion limit and the number of resolved 
sources were calculated for the Galactic disk binaries. 
However, for example. Galactic halo objects and extra- 
galactic binaries may also contribute to the GW signal in 
the LISA band. 
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The results of microlensing experiments can be con- 
sidered as evidence for the existence of massive com- 
pact halo objects (MACHO's). The most hkely MACHO 
mass i s between 0.15 and 0.9 M(T), depending on the halo 



keeps increasing at z ^ 3 (Kosenko & Postnov 1998). 
A computation with a star formation rate which is al- 
most constant at 1.5 ^ z < 5 and roughly the same in- 
put for stellar evolution as used in this study (Schneider 



model, and the total mass in MACHO's out to 50 kpc is 



9+;^ X 10 



^0: 



independent of the halo model (Alcock 



5+4 
et al 

(e.g. Gates et al. 1998), but two possibilities are relevant 
to this study. 



>000). The nature of MACHO's is still unknown 



The first is that they are white dwarfs ( Tamanaha et al 
1990|)~| A discussion is still going on whether the pres- 
ence of a significant white dwarf population is compati- 
ble with constraints derived from the chemical composi- 
tion of the halo and the cooling properties of white d warfs 
(e.g. IChabrier 1999| ; [Fields et al. 2000| ; [Hansen 2000|) . The 
fraction of binaries in this hypothetical population is un- 
known. However, because the star formation in the halo 



happened long ago (e.g. Adams fc Laughlin 1996 ), most 
short period double white dwarfs that could have formed 
will already have merged. For example, for a 1 Gyr long 
burst of star formation, 10 Gyr ago, with an IMF similar 
to the IMF of the disk, all halo close binary white dwarfs 
currently have orbital periods longer than ~ 0.3 hr, so they 
cannot contribute to the GW signal at frequencies higher 
than log/ ~ —2.75. The observed deficit of halo white 



dwarf progenitors in distant galaxies (Adams & Laughlin 
1996|)~^uggests that the IMF in the halo is peaked at or 
above 2 Mq , limiting a hypothetical double white dwarf 
population to even lower frequencies. Hence, we do not 
expect a change in the confusion limit due to halo white 
dwarfs, but they could contribute to the noise below this 
limit. 

Existing estimates of the co ntribution of halo do uble 
white dwarfs to the GW noise ( Hiscock et al. 200C ) are 
based on a simple rescaling to t he halo of the e stimate of 
the GW signal from the disk by Hils et al. (1990 ). Because 
of the different evolutionary histories of disk and halo this 
is unrealistic. Additionally, Hiscock et al. (200C ) use the 
lowest existing observational estimate of the local white 
dwarf space density, probably overestimating the relative 
importance of a possible halo white dwarf population. 

Another possibility is that MACHO's are low-mass 
black holes. Formation of low-mass (^ IMq) black 
holes is possible in inflationary cosmological models (e.g. 
Naselskii & Polnarev 1985). Further, as was shown by 
Nakamura et al. (1997), these black holes may form bi- 
naries. An estimate by Hiscock (1998) shows that under 
certain assumptions about the separations of the compo- 
nents, the GW background formed by halo binary black 
holes can be much stronger than the signal from the 
Galactic binary white dwarfs. If this model is correct, the 
noise produced by halo objects would bury virtually all re- 
solved signals from Gal actic systems (compare our Fig.^, 



et al. 200C), showed that only just above the point where 
the Galactic disk double white dwarf background drops 
sharply (around log/ = —2.75) the extragalactic back- 
ground could exceed the Galactic one. However, the sig- 
nals of the resolved binaries at these frequencies are at 
least an order of magnitude stronger than this background 
an d will probably be det ectable (see our Fig. |^ and Fig. 12 
of [Schneider et al. 200d|) . 

Finally it should be noted that, in addition to our es- 
timate above, a few tens of double neutron stars and neu- 
tron star - white dwarf binaries in globular clusters will 
probably be resolved (Bcnacquista et al. 2001). 



5.2. Comparison with previous studies 

The Galactic GW background produced by double white 



dwarfs was studied earlier by e.g. 


Evans et al. (1987); 


Lipunov et al. (1987 


); 


Hils et al. (1990) 


: 


Postnov & 


Prokhorov (1998); 


Webbink & Han (199^ 


i). The most 


widely quoted study is the one by 


Hils et al. (1990), who 



calculated the background based on the estimates of the 



number of systems by Webbink (1984 ). Because Webbink 



found a considerably higher birth rate of close double 
white dwarfs than we find. Hills et al. found a higher noise 



background. The same holds for the study by Evans et al 
|(1987| ), who in addition use a different Galactic distribu- 
tion. When we estimate the difference in the total number 
of systems in the Galaxy from the different birth rates, 
we find that the renormalised background levels differ by 
a factor ^ 3. 

Webbink & Han (1998) use a current birth rate of 
0.03 yr~^, similar to ours, but a constant star formation 
history and a larger age of the Galactic disk, of 15 Gyr 
( Han 1998 ). Their estimate of the background is slightly 
higher that what we find, probably due to the higher av- 
erage chirp mass (0.42) of their (wd, wd) systems or their 
assumed Galactic distribution, which is slightly different 
from our Eq. (^. Because they also calculated individ- 
ual systems (which they later average), they calculate the 
number of resolved systems (although with a different 
criterium). Above the resolution limit found by Han & 
Webbink (log / = —2.44) we have 3615 resolved systems 
in agreement with their number of 3600, despite differ- 
ences in the underlying white dwarf population. 

The fact that most studies give an estimate of the con- 
fusion limit within a factor ^ 2 from ~ 1.6 mHz found by 
us is a consequence of weak dependence of this limit on 



left panel and Fig. 2 of Hiscock 1 



the parameters of the models [see Eq. (19) of Evans et al 
This would resuh 1(19871 )]. 



in the non-detection of any resolved systems by LISA, and 
show up as an anisotropic noise. 

A significant contribution from extragalactic binaries 
to the background is expected only if the star formation 



The backgrou nd due to semi-detach ed white dwarfs 
was calculated by Hils fc Bcirder (200C| ). They conclude 
(as we do) that these stars are not important for the over- 
all background. This can be understood as a consequence 
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of their low strain amplitude due to their low chirp mass, 
which makes them unimportant even in the frequency 
range where they outnumber the detached systems. 

6. Conclusion 

We calculated the gravitational wave signal of Galactic 
disk binaries containing two compact objects. We discuss 
three populations: (i) double white dwarfs (including semi- 
detached systems) which produce a confusion limited noise 
background at low frequencies (log / ^ —2.8), (ii) resolved 
binaries and (iii) unresolved systems that have such a 
strong signal that they may be detected above the noise. 

The confusion limited background is dominated by de- 
tached double white dwarfs, although in a small frc^qucncy 
range (—3.4 ^ log/ ^ —3.0) the semi-detached systems 
form the majority of the systems in each bin. The dou- 
ble white dwarf gravitational wave background, which in 
our model consists of the sum of the signal of all 150 
million systems in the Galaxy shows large spikes caused 
by strong-signal (i.e. close) systems, which might be de- 
tectable above the noise. 

Adding the binaries containing neutron stars and black 
holes (which are much less numerous than white dwarf 
pairs), we find the distribution of bins containing one, 
none and more than one system and show that the 
"confusion limit" as a single value does not exist: at 
log / = —2.84 the first resolved bins are found, while up 
to log / = — 2.3 bins containing more than one system are 
present. 

We find 12171 resolved systems of which the vast ma- 
jority are double white dwarfs. There are only 8 double 
neutron stars and 38 neutron star - white dwarf binaries 
resolved. Finally we calculate that there are 6104 systems 
(5943 double white dwarfs, 124 neutron star - white dwarf 
systems, 31 double neutron stars and 6 systems containing 
a black hole) which have a signal well above the double 
white dwarf background and the LISA sensitivity level and 
might be detectable. 

Out of 12124 resolved double white dwarfs, 501 have a 
combined mass above the Chandrasekhar limit and peri- 
ods short enough to merge in 10 Gyr and are thus potential 
SN la progenitors. Such double white dwarfs have not yet 
been found optically. If a system would chirp, LISA will 
measure the chirp mass and the distance to the system 
allowing a good estimate of its total mass. The number 
of chirping SN la progenitors is ~ 94 for an integration 
time of 1 yr. But even for these systems the actual time 
to coalescence is ^ (4500 - 80000) yr. 
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Appendix A: The formation of black holes 

The simplified description of the evolution of massive stars 



used in thi s paper is based on results of Eggleton et al 
1989D and [Schaller et al. (1992|) and can be summarised 



as follows (see also Portegies Zwart & Yungelson 1998) 



1. The radii of the massive stars are limited to 1000 R 





(after [Schaller et al. 1992| ) . 

The maximum amount of mass loss for hydrogen rich 
massive stars is given by 0.01 M^, unless the whole 
envelope is lost and the star becomes a Wolf-Rayet 
star. Thus stars above 85 Mq evolve off the main se- 
quence immediately to become Wolf-Rayet stars with 



an initial mass of 43 Mq (after Schaller et al. 1992). 
Wolf-Rayet stars lose mass according to the equa- 
tion proposed by [Langer (19"89| ): M = 5 x IQ-'^M'^-^ 
Mq yr-1. 

In the supernova 50% of the m ass of the exploding ob - 
ject is ejected (consistent with Nelemans et al. 1999 ). 
Exploding helium stars more massive than 10 Mq col- 
lapse to a black hole. Thus, the lower limit for black 
hole masses of M = 5Mf7^ is consi stent with observa- 
tional estimates (e.g. Charles 1998[ ). 
Black holes do not get a kick at birth in cont rast to 
neutron stars (e.g. White & van Paradijs 1996 ). 



